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Abstract 

We discuss several new results on nonnegative approximate control- 
lability for the one-dimensional Heat equation governed by either multi- 
plicative or nonnegative additive control, acting within a proper subset of 
the space domain at every moment of time. Our methods allow us to link 
these two types of controls to some extend. The main results include ap- 
proximate controllability properties both for the static and mobile control 
supports. 
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1 Problem setting 

In this paper we study the approximate controllability properties of the one 
dimensional Heat equation governed either by an additive nonnegative locally 
distributed control or by a multiplicative control of local support. For simplicity, 
we will assume that the space domain is the interval (0, 1). More precisely, we 
consider the following two Dirichlet problems: 
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{yt{x,t) = y^^{x,t) +u{x,t)Xoj{t){x), in Qt = (0, 1) x (0, T), 
y{Q,t) = y{l,t)=Q, in (0,r), (SI) 

y{x,0) = yo{x), in (0,1), 



and 

yt{x,t) = y^^{x,t) +v{x,t)y{x,t)xuj(t)ix), in Qt, 

y(0,t) = yil,t) = 0, in (0,T), (52) 

y{x,0) = yo{x), in (0,1), 

where yo e i+(0,l) = £^(0,1) : </>(.x) > a.e. in (0,1)}, with yo ^ 

for {S2). In the above, u is the additive locahy distributed control, while v is 
the multiplicative (or bilinear) control (in their respective systems), each one 
supported at each time f in a proper subinterval uj{t) of [0, 1]. We will assume 
that 

cj{t) = ir{t),r{t)+l), (1.1) 

where the length I e (0,1) is fixed and r : [0,T] — > [0, 1 — Z] is a function 

that defines the position of the support at time t. We call the support oj{t) 
static if it does not change in time and mobile if it does. The symbol Xw(t) 
denotes the characteristic function of u){t) (i.e. Xu{t){x) = 1, when x G u){t) and 
Xui{t){x) = 0, for X ^ w(t)). In the sequel, we will denote by fjy^^u.r the solution 
of system (5*1) and by yy„,v,r the solution of system (5*2). By the linearity of 
(5*1) it is clear that yyg,u,r = jjyo +yo,u,r, where yy^ denotes the unique solution 
of (5*1) with u = 0. If the subinterval uj{t) = uj for all t, (i.e. it is static), we 
will use the notation yyg^u,cu instead of yyo,u.r (rcsp., yy^^v.u instead of yy^.v^r)- 

The control u in (5*1) represents an external source acting within uj{t). In 
turn, the multiplicative control v in (52) can be viewed as controlling the reac- 
tion rate of the process (52) within its support. 

We will investigate how rich the sets of solutions to the above controlled 
systems are at the fixed final time T > 0. To that end, let us introduce the 
following notations: 

Fi,yo,mbiT) = {yy„uAx,T) : X e (0, 1), u e LHQt) and r e PC[0,T]}, 



F2,yo,mb{T) = {yyo,vAx,T) : X £ {0, 1), V e L°°(Qt) and r e PC[0,T]}, 

(1.2) 

where PC[0,T] denotes the set of functions r : [0,T] — > [0,1 — I] such that 
r is picccwisc constant with (at most) a finite number of discontinuities and 
LHQt) ^{(f>e L\Qt) : 0(x, t)>0 a.e. in Qt}. 

System (51) is said to be approximately nonnegative controllable at time T > 
if for any given ya G Ll(0. 1), the set F+y^^^^^{T) is dense in ?^(-, T)+Ll{0, 1). 
Taking into account the linearity of the PDE, it is quite easy to verify 

Kyo,mbiT) = yy,{., T) + F+^^^.iT). (1.3) 
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Hence, the approximate nonnegative controllability at time T > is clearly 
equivalent to the density of F^^ mfc(^) -^+(0^ !)• Analogously, system {S2) is 
said to be approximately nonnegative controllable at time T > if for any given 
yo e i+(0, 1), yo ^ 0, the set F2,y„,mb{T) is dense in L+(0, 1). For yo = 0, by the 
uniqueness of solution for the system (<S'2), it is obvious that F2fi^mb{T) = {0} 
and the aforementioned density property in i+(0, 1) clearly does not hold. 

Formally, it is apparent that the PDF in (32) can be obtained from the 
corresponding one in {SI) by setting 

vix,t) = ^ (1.4) 
y{x,t) 

This identity will be a key point to get our results. Of course, these two problems 
are not equivalent in the sense that in (51) the control u is always nonnegative, 
while V in {S2) can be negative as well, thus providing a "qualitatively richer" 
set of solutions for (52) compared to that of (51). 

When the support of the control is the whole space domain, it was shown in 
[8] that the set 

Fi\,iT) ^ {yy^,^^Jx.T) : x E (0,1), u e LHQt), c^=(0,l), Vte(0,T)} 

(1.5) 

is dense in yyg{-,T) + L^{0, 1). Furthermore, it was shown in [9]-[ll] that the 
corresponding property also holds for the system (52): i.e., the set 

F2,yjT)^{yy,,,,Ux,T): xeiO,l), veL°-iQT), ^=(0,1), Vie(0,T)} 

(1.6) 

is dense in L^{0, 1). In fact, these results remain true for more general linear 
and semilinear parabolic equations in several space dimensions. In this paper, 
however, we are specifically interested in the case when uj(t) is a proper subin- 
terval of (0, 1) at every moment of time. 

Let us emphasize that the multiplicative controllability problem is of non- 
linear nature, even in the case of the linear equation like (52), because the 
mapping v — ^ yyo,v.u> is highly nonlinear: for instance, when w = (0, 1) and v is 
constant, it is well known that yy„.v,u]{x,t) = yyi^(x,t)e~^* . This means that we 
cannot use here the traditional duality argument, which is the typical strategy 
to approach controllability problems in the case of linear PDF with additive 
controls. 

Among earlier works on the controllability of linear PDF by means of bilinear 
controls we refer to the pioneering work [3j by Ball, Mardsen, and Slcmrod. The 
paper [13] further explored the ideas of [3]. 

An extensive and thorough bibliography on controllability of bilinear ODF 
is available. Let us just mention the survey [1]. Research in this area was 
seemingly originated in the 60s, on the one hand, by the works of Kucera, 
who linked this area to the Lie Algebra approach, and, on the other hand, by 
the works of Rink and Mohler [20] (see also [12]), who pursued the qualitative 
approach and numerous applications. 
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For recent works on global controllability of parabolic and hyperbolic linear 
and semilinear PDE we refer to [6], [8]- [10] and specially to the monograph [Tl] 
and the references therein. 

We also refer to the works (see also the references therein) on the 

issue of optimal bilinear control for various PDE. A closely related issue is sta- 
bilization by means of bilinear controls, see [2] and [TB] • 

2 Main results 

First of all, we will show that for static subdomain the above controllability 
property is out of the question for both systems (SI) and (5*2). 

THEOREM 1 (Lack of controllability for static support) Let us assume 
that Ho G ^+(0, 1), j/o ^ 0, and uj is a fixed proper subinterval of (0, 1). Then, 
for any T > 0, the sets 



F+o^{T) = {yo,uA^,T): xe{0,l), ueLliQT), u{t) = u:, Vte(0,T)}, 

(2.1) 



F2^y,AT)^{yy,,vAx.T): x€{Q,i), V ^l°°{Qt), uj{t) = uj, vte(o,r)} 



are not dense in i^(0, 1). 

Nevertheless, we can prove (in a constructive way) the following ^HocaF 
controllability result. 

THEOREM 2 (Local nonnegative controllability) Let us assume that ui 
is a fixed proper subinterval o/ (0, 1). Then, for any T > 0, the set 



Pi,ojT)\^ = {yo,uAx, T): x£uj, ue L^Qt), io{t) = uj, Vi G (0, T)}, 



is dense in L^(w). 

If the control support is allowed to move with an additional piecewise con- 
stant control r{t), defining its position at time t, we have the following global 
controllability result. 

THEOREM 3 (Approximate controllability with mobile support) Let 

us assume that yo G i+(0, 1), yo ^ 0. Then, for any T > 0, the sets F^^ mbC^) 
and F2^y„^mb(T) are dense in L^(0, 1). 

In the proof of Theorem [T] we have faced with the following system: 



and 



(2.2) 



(2.3) 



yt{x,t) 



{x,t), in Qt, 

uo{t), in (0,T), 

u,{t), in (0,T), 

yo{x), in (0,1), 



y(i,i) 

. 2/(x,0) 



{S3) 
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with uq and ui acting as boundary controls, meanwhile i/q is fixed. 
If we introduce the set 



FyoAT) = {y, 



{x,T) : xe{0,l), uo,ui€L^{0,T)} 



(2.4) 



where yyo,uo,ui denotes the unique solution of (S'3) defined by transposition (see 
[171 pp. 194-7]), we have derived a related result on the lack of nonnegative 
approximate boundary controllability with nonnegative boundary controls, in 
the following terms: 

THEOREM 4 Let us assume that ya e L+(0, 1). Then, for any T > 0, the 
set Fyg,i,{T) is not dense in yyg{-,T) + L+(0, 1). 

An essential tool for obtaining these results are the classical Maximum Prin- 
ciples for parabolic equations (both in the weak and strong form). Since they 
will be cited very frequently along the paper and there are different versions 
depending on the framework, we recall them here: 

THEOREM 5 (Maximum Principles) Let us assume that yo e L^(0,1), 
u € L^{Qt) and v G L°°{Qt). Let y G L^{0, T; H\0, 1)) n C([0, T];L^{0, 1)) be 
a solution of the problem 

( ytix,t) = yxxix,t) +v{x,t)y{x,t) +u{x,t), in Qt, 



i) If yo > in (0,1), u> in Qt, y{0,t) and y{l,t) > in {0,T), then 
y{x,t) > a.e. in Qt- 

a) lfu=0,v<0 in Qt, y{0, ■),y{l, ■) e L°°{0,T) and yo G L°°{0,1), then 
y{x,t) < K a.e. in Qt, with 

K = max{||(j/(0, •))^IIl-(o,t), ||(2/(1, ■))^IIl-(o,t), II (2/o)^||l~(o,i)}' 

where z+ — max{z,Q\. 

m) (-Strong; //ye 2.2(0, T; iJi(0, 1)), yo > Q m (0,1) (yo ^Q),u>Qin 
Qt, V = v{x) <0 in (0, 1), then y{x, t) > in (0, 1) for all t € (0, T] . 

For the proofs, we refer (among others) to [Ml Chapter III, section 7, pp. 
181-191] and W, pp. 21-23]. See also T^, Chapter I, section 2, pp. 11-25] for 
related results in the classical framework. 



We begin with the proof of Theorem[T]for system (51). Our goal is to construct 
a simple explicit counterexample. 




in (0,1). 



3 Proof of Theorems [T], H and [H 
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3.1 Proof of Theorem [T] for system (SI). 

Since u g L'^{Qt), it is well known (see for instance [HI Theorem 4.2, p. 
160]), that each solution ?/o,„,^ belongs to L^{0, T; H^{0, l))nC([0, T];L^{0, 1)). 
Moreover, due to the Theorem [5]-i) 

yQ,u,u{x,t)>0 a.e. x£ (0,1), Vie [0,T]. 

Consequently, F^^f^ ^{T) C i+(0, 1). We will show that F^^(r) is not dense 
inL^(0,l). 

Since w is a fixed proper subinterval of (0, 1), there will exist a sufficiently 
large natural number m > 2 such that ut C (^,1). Let us fix that m for the 
rest of the proof. It is now straightforward to check that if 

p{x,t) = ^{x)exp („,_i)2 
with 

J sin(m7ra:), a; G [0, :^), 

^^'''-X (l_^)sin(^(l2^), xe[i,l], 

p G C^{Qt) and it is the unique classical solution of the adjoint problem 

-pt{x,t) = pxxix,t) + h{x,t), IuQt, 
piO,t) = p{l,t) = 0, in (0,r), (3.1) 

p{x,T) = ifi{x), in (0,1), 



where 



{(m-2)m^7r^ • / \ /m^7r^(t-T)\ ^ rf, 1\ 

sm (mTrx) exp ' j, x€[0,-), 

0, ^e[i,i]. 



Furthermore, 



h{x,t) > in Qt, (3.2) 
1 



,1 
m 



x[0,r]. (3.3) 



p{x, t) < in 

Let us now pick up the state in L^{0, 1) 

yd{x) = max{^(a;),0}. 

We will see that yd is "unachievable". Arguing by contradiction, let us 
suppose that there exists a sequence C L^{Qt) such that: 

yo.Uk,Lj{-iT) — > yd inL^(0,l) as fc ^ +oo. 

Multiplying (51) (with u — Uk, i^{t) — w and j/o — 0) by p, further integra- 
tion by parts, combined with p.l^ - p.Sp provides: 
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yo,uk,ujix,T)ip{x)dx < / yo^uk,coix,T)ip{x)dx+ I h{x,t)yo^uk,uj{x,t)dxdt 



10 Ja 



Ukix,t)p{x,t)dxdt < Vfc. (3.4) 

uix{0,T) 

Thus, we arrive at the contradiction 

0< / {s'm{Tmrx))'^dx — / yd{x)(p{x)dx — \im / yo,uf.^u:{x,T)ip{x)dx < 0. 
Jo Jo ^ Jo ' 

This ends the proof of Theorem [T] for system (5'1).B 
3.2 Proof of Theorem [21 

It is enough to show that solutions to system {SI) (with j/o = and fixed w) at 
time T can approximate in L^(a;) any element of Hq(uj) P| L^(a;). This can be 
shown by a simple modification of the proof given in [S] for the case lj — (0,1). 
Nonetheless, the argument used there does not provide any insight on what 
controls are used to achieve the desirable steering. Therefore, here we give a 
different constructive proof of Theorem [21 through Fourier series expansions, 
which can be used to prove the aforementioned result in [5] as well and will also 
serve us as a tool to prove the Theorem [3| below. 

Consider any element yd G ^^0(1^) Pi Extend it by zero to the whole 

interval (0, 1) and, for simplicity of notations, denote the resulting function by 
yd again. Clearly, it lies in -ffg (0, 1) p| i^(0, 1). Given any 5 € (0, T), we define 

0, te[o,r-5), 

1 

<5 

Taking now advantage of the Fourier series expansion, we know 



Vd 



(x) = 2^afesin(fc7rx), inHl{0,l), (3.6) 



k=l 



with Ofe = yd{s)sm.{k-Ks)ds. 

The classical method of separation of variables provides us the following 
expression for the corresponding solution of (5*1) with yo = 



0, te[o,r-J], 

' '-"'^l^r^'^ akSinikTrx), t e [T - S,T]. 

Hence, we derive that 

11 yo.usA;T)-yd 11^1(0,1) -2^' E«fc^' ( ^2fc2j 1) ' (3-8) 



7 



where we endowed the Sobolev space Hq{0, 1) with the norm 

1/2 



Let us show that the right-hand side of p.Sp tends to zero as (5 — > 0+. To 
that end, we introduce the auxihary function 

V'(r) = (^^—^ 1^ > r>Q. 

It can be verified in a straightforward way that ?/; is a C°° function, strictly 
increasing in (0, +00), ?/'(0) = and V'([Oj +00)) C [0, 1). 

Taking into account that yd G 7?q(0,1), given any e > 0, we can find a 
natural number such that, 

00 

E -Ik' < (3-9) 

k=N+l 

Using now that ip is continuous at and "0(0) ~ 0, we can select j > small 
enough to guarantee that 

N 

^{7T'N^)Y,4k' < V<5e[0,5]. (3.10) 

fc=i 

Combining previous estimates with the fact that ?/j(7r^fc^(5) < 1 for all k, we 
get that 

CN 00 \ 

k=l k=N+l I 

(3.11) 

for all b e (0,(5]. This ends the proof of Theorem [21 but for proving later 
Theorem [3] corresponding to system (5*1), it is convenient to make a small 
modification to previous construction as follows: instead of the control us{x,t) 
given by (|3.5|) . we can take a control with value zero at the end of the process. 
More precisely, we can consider for instance 

0, te{0,T-25), 
usix,t) = { ^dix), te[T-2S,T-5), (3.12) 
0, te[T~5,T]. 

Applying the method of separation of variables in each time interval, it is 
not difficult to see here that the solution of (SI) with j/g = is given by 
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r 0, 



k=l 



Ofc sin(A:7ra;), 



t e [0,T-2S], 

te [r-2(5,r-(5], 



'^k^S 2,2/ 



2 v^oo 1-e 

and therefore, arguing as before, we arrive to 



fc (*"r+'5)afcsin(fc7ra;), t e [T - S,T], 

(3.13) 



2/0, «5 



fc=i 



(3.14) 

We can finish as above, seeing that the right-hand side of p.l4p tends to 
zero as (5 — !■ 0"*" with the help of the new auxihary function 

^-(0 = 1 , ^ > 0, 



that satisfies the same properties than tpM 

3.3 Proof of Theorem [1] for system (5*2). 

Let (a,/3) be any subinterval of (0, l)\w and let us pick out any yd G -^+(0, 1) 
such that Hd = in {a,P). Arguing by contradiction, let us assume that there 
exist {wfc}fcLi C L°°{Qt) such that 

yyo,vkA'^T) ^ Vd ini+(0, 1) asfc^+oo. 

In particular, we have 

2/j/o,«fc,c^(-,?')|(„,^) ^ inL^(a,/3) as fc ^ +oo. 

The restriction of each yy„ ^vk.uj to the interval (a,/3) can be viewed as the 
solution of one system of the type 



yt{x,t) = yxx{x,t), in (a,/3) X (0,T), 

y{a,t) = UQ{t), in (0,r), 

y{l3,t) = u,{t), in (0,r), 

y{x,0) = yo{x), in {a, (3), 



{S3a) 



with nonnegative boundary values at the extremes a and /3, because yyQ^v^,u is 
nonnegative thanks to the Theorem [5]-i). 

On the other hand, if we denote by ij the unique solution of the system 



yt{x,t) = yxx{x,t), in (a,/3) X (0,T), 
y{a,t) = y(/3,t) =0, in (0,r), 
y{x,0) = yo{x), in {a, (3), 



(536) 
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as a consequence of the Theorem O-i) we have that 

- y{x,t) > 0, a.e. {x,t) G {a, /S) x (0,T), 
and thanks to Theorem [5]-iii) 

y{x,T) > 0, a.e. x G (q;,/3). 
The contradiction arrives combining ah the previous conditions in the form 
= lim \\yy^,v^A-^T^)\\h^{a3) ^ lli/(->^)llL2(a,,3) > O.B 

3.4 Proof of Theorem HI 

In the previous proof, we have faced with a problem of the type (S'3) (more pre- 
cisely, (5*30)). Now, we will show that the nonnegative controllability property 
for these problems is out of the question, even if the obstruction phenomenon 
due to the initial datum is removed. 

By the linearity of the problem, it is well known that yyg^uo.ui — Vyo +yo.ua,ui , 
and the following relation holds 

Fy,AT)^yyo{-,T)+FoAT). (3.15) 

Therefore, the conclusion of Theorem S] is equivalent to say that Fq i,{T) is not 
dense in i+(0, 1) and this is what we will prove. 

At this point, let us remind that (in general) the solution yo,uo,Mi niust be 
defined by the transposition method (see [171 PP- 194-7]). In our situation, its 
expression is given by 

00 

yo,uo,udx^'t) =Ybk{t)V2sm{kTTx), (3.16) 

/c=l 

where 

bkit) = V2kTT e-^'"'(*-^)(wo(s) + (-l)^'+'Mi(s))ds. (3.17) 
Jo 

It is not obvious at all that j/o,tio,«i ^) ^ L^{0, 1)- In fact, in [T71 p. 202], it is 
constructed an example where yo,uo,«i 

(•,r) ^ L^{0, 1), by taking uo(t) = and 
ui{t) = The difficulty here is that ui e L'^{0,T), but not in L°°(0,r). 

In our case, using (I3.16p - (I3.17|) . we can estimate 

/ l2;o,«o.«i(a:,T)l2da; = ^15fc(r)l2<CiV-2 <c», (3.18) 

fc=i 

with Ci = (I|mo1Il~(o,t) + I1wi1|loo(o,t))^- 

Moreover, applying the Theorem[5]-i) we derive that yo.uo,ui{x,T) > 0, a. e. 
X € (0, 1). Consequently, Fo^f,(T) is included in i^(0, 1). Let us show that it is 
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not dense in i^(0, 1), arguing by contradiction. To this end, it can be checked 
by direct calculation that 

p{x, t) = - exp (97r2(t - T)) sin (Sttx) 

is the unique classical solution of the adjoint problem 

-Pt{x,t) = p^^{x,t), in Qt, 

p{0,t) = p(l,t) = 0, in (0,T), (3.19) 
p(x,T) = — sin(37rx), in (0,1). 

Moreover, it satisfies p2;(0,i) < 0, Px{y,t) > in (0,T). Now, we select 

2/d(a;) = niax{p(x,r),0}, a;e(0, 1), 

as the target state in L+(0, 1). 

If Fo^b{T) is dense in i^(0, 1), there exist two sequences {uofc}fcLi, {""ifcjfcLi 
in L'^{0,T) such that for the corresponding solutions {Uk} kLi of (53) with 
2/0 = (namely, yk = 2/o,«ofc,"ifc): '^e have: 

2/fe(-,T) — >yd in L^(0,1) as fc ^ +0O. 
Multiplying by p the PDE satisfied by yk and integrating by parts, we get: 

yk{x,T)p{x,T)dx = uok{t)Px{0,t)dt - uik{t)px{l,t)dt < Q. 



(Alternatively, previous equality can be derived from the expressions (|3.16p - 
p.l7p for each yk-) Thus, we arrive to the contradiction 

2/3 ,1 ,1 

(sin (37rx))^(ia; = / yd{x)p{x,T)dx = lim / yk{x, T)p{x, T)dx < O.M 



1/3 Jo k^°°Jo 

4 Proof of Theorem [3]. 

4.1 Proof of Theorem [3] for system (SI). 

To show that F^^ mbC^) dense in i+(0, 1), we will use the technique developed 
in the proof of Theorem [2] First, let us fix M the smallest natural number 
satisfying M ■ I > 1, where I is the (fixed) length of the mobile subinterval uj{t) 
(see (dH)). 

Given any element yd G L^(0, 1), we can decompose it into M "pieces" 
localized in disjoint subintervals of length (at most) I, as follows: 



Vdi (x) 
yd, (x) = 



yd{x), X e (0,0, 
0, xe{l,i), 



Vd{x), X e {{j - 

0, otherwise. 
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for j = 2,...,M-l, 



M - i y^'^'^^' xem ~ 1)1,1), 

ydM(x) - \ 0^ otherwise. ^ ^ 

Clearly, each Udj & L^(0, 1) and 

yd{x) = ^ydj{x), a.e. xe(0,l). 

Moreover, given e > 0, there exist ijdj € i?o(0, l)nL^(0, 1) for j = 1, . . . , M, 
with its support ujj contained in the corresponding one of each ydj and verifying 

M 

II i/d, - yd, ||l2(o,i)< 2' ("^-^^ 

Now, we can argue "by parts", using the constructions p.l2p and p.5p 
established in the proof of Theorem [2l as follows: given yd^ , there exists Si S 
(0, ^) such that if we take the control 

r 0, t G (0,T-2(5i), 

ui(x,t) = <^ (r-2,5i,T-<5i), (4.3) 

[ 0, t e (T-5i,T), 

it is satisfied 

II 2;o,«i,a.i(-,7') - ydi ||l2(o,i) < — . (4.4) 

Arguing now with y^j, there exists 62 G (O, ^) such that if we take the control 

r 0, teiO,T- 262), 

U2{x,t) = \ j^ydAx), te{T~2d2,T-62), (4.5) 
[ 0, teiT-d2,T), 

it is satisfied ^ 

II yo,u2,'^2i-,T)-yd2 \\l^o,i) < (4.6) 

Iterating this process, we obtain 6j e ^0, , j = 1, . . . , M, and controls 

r 0, (0,T-2<5,), 
^.(^,0 = < iy<i.(^). te{T-2S„T-Sj), (4.7) 
I 0, (T-(5„T), 



j = 1, . . . , M - 1, and 



0, te{Q,T-SM), 

<5a 



UM{x,t)^i ^^...(.), ie(T-5M,r), (4.8) 



for which 

2M 



2/o,«,.^,(-,r) -y<jJU2(04) < j = l,...,M. (4.9) 
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Wc finish the proof by taking 



M 



(4.10) 



and noticing that by the hnearity 



M 



Af 



y(- 



{■,T)-yd \\l^{o,i)< II \\lHo,i) + II ^'^^ ^ ^/d IIl2(o,i)< 



< X! II ^o,fy,a;j(-,T) - yd, ||l2(o,i) +^ II yd, ||l2(o,i)< e. (4.11) 

Let us remark that in each time subinterval {T — Sj,T — Sj^i), the value of 
the control u{x,t) coincides with the value of Uj+i{x,t) and hence its support 
is included in uj{t) = {r{t),r{t) + 0, with r{t) G PC[0,T] given by 



r 0, 

I, 
21, 



r{t) = <^ 



t € {T - 5i,T ^ 82), 



(4.12) 



iM-2)l, te{T-SM-2,T-6M-i), 
^ l-l, te{T~SM-uT). 

This implies that u{x,t) — u{x,t)xuj(t){x) and y = yo,u,r, as desired.l 



4.2 Proof of Theorem [3] for system (5*2). 

We will prove the result in two steps. First, following the notations of Section 
1 (compare with (|1.2p ). we introduce the set 

F2,yo,mbiT) = {yyoM^,T) I X e {0, 1), V G (Qt) and r e PC[0,T]}. 

(4.13) 

Applying twice the Theorem [5]-i) (first to yy„,v,r and later to the difference 
yyo,v,r — yyoj treating the bilinear term as a free one), it is well known that 

yyo,vAx,T) > yy„{x,T), a.e. x £ (0,1), Vv e L°^{Qt), Vr e PC[0,T]. 



4.2.1 Step 1. 

We will prove that the system (52) is approximately nonnegative controllable 
at time T, in the following sense: 



13 



THEOREM 6 (Nonnegative multiplicative mobile control) Let us as- 
sume that 2/0 G ^+(0, 1), 2/0 ^ 0. Then, for any T > 0, the set F^y^^ mbC^) 
dense in yy„{-, T) + i^(0, 1). 

Proof of Theorem [6l Pick out a generic element 2/d € ^+(0, 1). As a conse- 
quence of Theorem [3] for system (5*1), given e > 0, there exist u £ L^{Qt) and 
r € PC[0,T] such that 

II yQ,u,r{-,T) - yd ||l2(o,i)< e- (4-14) 

By weU known density results, we can suppose, without loss of generality, 
that u G C-^.{Qt) (that is, u is continuous and nonnegative in Qt) with compact 
support K included in (0, 1) x {S,T) C Qt, for some 6 > 0. 

Using classical regularity results (see [71 Theorem 10.2, p. 140] and [l4l 
Theorem 9.1, p. 341-2] it is also known that jjy^ € C°°(Qt) n C+ ([0, 1] x [S,T]). 

Now, let us consider y — yy^.u.r and 

v(x,t)^i l^^''^^^^' (4.15) 

\ 0, {x,t) eQT\K. ' 

By the Theorem [S]-i), it is clear that y > 0, a.e. in Qx- We will prove 
below that v G L°^{Qt)- Using this fact, it is straightforward to see that 
y = yyo,v.r (i-e. y is the solution of system (5*2) corresponding to 2/o,w,r), 
because u{x,t)xuit){x) = v{x,t)y{x,t)xuj{t){x) in Qt, and 



II yyo.v,r{-,T) - {yyo{-,T) + yd) Wmo.i}^ 

= 11 yyo,n,r{-,T)~{yyg{-,T) + yd) ||L2(oa) = || yo,u,r{-,T)-yd \\moA)< e, (4.16) 

as we were looking for. 

Hence, it remains to show that v £ L°^{Qt)- Since jjyg e C+([0, 1] x [5,T]), 
Vyo > in Qt and X is a compact set, there is a constant p > such that 
Vyo > p in i^. In fact, we can guarantee that p > 0: if there exists {xi,ti) € 
K C (0, 1) X {S,T) such that yyg{xi,ti) = p = 0, thanks to the Theorem [5]-iii), 
we conclude that jjy^ = in [0, 1] x [0,ti] and, consequently, 2/0 = in (0, 1), 
which it is excluded by hypothesis. Then, we get 

II II , 11 " IU°°(g) _ 

II V ||l~(Qt)< •■ 



4.2.2 Step 2. 

We want to approximate any element yd e L^(0,1), even if yd ^ yyg{-,T) + 
L^(0, 1). Roughly speaking, the main idea here is to take very large negative 
controls at the beginning of the process and move the control support across 
the space domain, in order to reduce the lower barrier that represents yyg{-,T), 
making it small enough. 
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For that purpose, we fix M (as previously) the smallest natural number 
satisfying M ■1>1, where I is the length of the mobile subinterval uj{t). Given 
some intermediate times < Ti < T2 < . . . < Tm < T, let us denote by yi{x, t) 
the unique solution of the following system 

yt{x,t) = y^^{x,t) -■miy{x,t)x{o,i){x), in (0, 1) x (0, Ti), 
y(0,t) = y(l,t)-0, in (0,Ti), (4.17) 

y{x,0) = yo{x), in (0,1), 

where mi > will be carefully chosen later. Assuming now that yj^i is already 
known, yj will denote the solution of the system 

{yt(x,t) = ya:x{x,t) ~mjy{x,t)x{(j-i)iji){x), in (0, 1) x (Tj^i, T,), 
y{0,t) = y{l,t)=0, in (r,-i,T,), 

y(x,Tj_i) = yj^i{x,Tj^i), in (0,1), 

(4.18) 

with nij > for j = 1, . . . , M — 1 also to be chosen. Finally, yM will denote the 
solution of the problem 



yt{x,t) = yxx{x,t) ~ mMy{x,t)x((M~i)i,i){x), y. {Tm-i,Tm), 

y{Q,t) = 2/(1, t) =0, in {Tm-i,Tm), 

y{x,TM-i) = yM-i{x,TM-i), in (0,1), 

(4.19) 

with mjv/ > to be chosen. The following result will be essential in this section: 

THEOREM 7 With previous notations, given e > 0, there exist nonnegative 
constants mi, . . . ,mM and intermediate times < Ti < T2 < . . . < Tm < T 
such that ^ 

|lyM(-,TM)||L2(o,i) < -. 



We postpone its proof to the next sections, but we will conclude now the proof 
of Theorem [21 Given yi{x,t) the unique solution of the system 

yt{x,t) = yxx{x,t), in (0, 1) X (Tm,T), 
yiO,t) = y(l,t) = 0, in iTM,T), (4.20) 
y{x,TM) = yM{x,TM), in (0,1), 

it is well known that 

l|yi(-,r)|U.(o,i) < ||2/M(-,rM)||L2(o,i) < |, (4.21) 

where the last inequality comes from Theorem [71 

Thanks to the Theorem [5]-i) , it holds that yj{x,t) > a.e. {x,t) S Qt for 
j — 1, . . . , M. In particular, ynix, Tm) > a.e. x G (0, 1). 

Assume that yd € L^(0,1) and e > are fixed. Applying now Theorem [6] 
for the corresponding system of type (52) in the domain (0, 1) x {Tm,T) with 
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initial datum yM{-,TM), we deduce the existence of V2 & i+ ((0, 1) x {Tm,T)) 
and r2 G PC[Tm ,T] such that the solution y2{x,t) of the system 



yt{x,t) = y^^{x,t) +V2{x,t)y{x,t)xu(t){x), in (0, 1) x (Tm, T), 
2/(0, = 2/(1, = 0, in (rM,T), 

y{x,TM) = yM{x,TM), in (0,1), 



satisfies 



II y2(-,T) - {yd + Vi{-,T)) ||l2(o,i)< -. 



We finish the proof, by noticing that 

yyo,v,r{x, t) - 



with 



yi{x,t), e (0,1) X (0,Tm], 

y2{x,t), {x,t)e (0,1) x[Tm,T), 



-mi, (a;,i) e (0, 1) X (0,Ti), 
-TO2, (a;,t) e (0,1) X (ri,T2), 



v{x,t) = { ~mj, (a;,t) e (0,1) X (r,-_i,Tj), 



-TOM, {x,t) e (0,1) X (rM-i,rM), 

[ V2{x,t), {x,t) e (0,1) X {Tm,T), 



and 



r(t) 



r 0, 



te{Ti,T2), 
t e (TjjTj+i), 



(M-2)?, t e (TAf-2,rM-i) 

1 — ^ i ^ {Tm-i,Tm), 

[ r2{t), te{TM,T). 



(4.22) 
(4.23) 

(4.24) 



(4.25) 



(4.26) 



Clearly, v £ L°°{Qt) and r e FC[0,T]. Moreover, thanks to (|4:2T1) and 
(j4?23l) . it is satisfied 

II yvo,v,r{-,T) - yd ||l2(o,i) = || yya^vA--, ^) " (j/d + J/lC'' ^)) + J/lC': ^) IU2(0,1)< 
<|| yyo,v^r{-,T) - (yd +yi(-,r)) ||l2(0,1) +I|2/i(->T')IIl2(04) = 

= 11 y2(-,T)- (2/d + yi(.,T)) |U.(o,i) +l|yi(-,7^)llL^(04) < (4-27) 
as we were looking for.B 
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4.2.3 Auxiliary results. 

Our proof for Theorem [7] is quite technical and relies on the following results: 

PROPOSITION 1 Assume that v e C[0, 1] with v{x) < for all x E [0, 1] 

and yo G C^+^p, 1], with a S (0, 1) verifies yo{0) = yo{l) = y^,'(0) = = 

and yo{x) > for all x S [0, 1]. Then, the classical solution of the problem 

ytix,t) yxx{x,t) +v{x)y{x,t), in Qt, 

y{0,t) = 2/(1,0-0, in (0,r), (4.28) 

y{x,0) = yo{x), in (0,1), 

satisfies the following properties: 

a) yt{x,t) < maXxeio,i]{ya{x)y, V {x,t) G Qt, 

b) 0<2/,(0,i) <e-max,g[o,i] (2/0(2^)6"°^"^), Vie [0,r] and 

c) e-min,e[o4](yo(a^)e^''^"^)<y-(l'0<0, Vte[0,r]. 

Proof of Proposition [H In the case yo = the result is clearly valid, 
because y = 0. So, let us suppose yo ^ 0. 

The existence and uniqueness of a classical solution y e C^+'^'^+'^/^(Qt) for 
problem (|4?28| is a consequence of [14( Theorem 5.2, p. 320]. Moreover, by the 
Theorem [5]-i) and ii), it is also known that 

< y{x,t) < ||?;o||l~(o,i) ^ (^^'0 G Qt- 

To prove a), we differentiate (|4.28p with respect to the time variable to 
deduce that yt is the solution of the problem 

zt{x,t) = Zxx{x,t) + v{x)z{x,t), in Qt, 
ziO,t) ^ z(l,t)=0, in (0,r), (4.29) 

z{x,0) = y'oix) +v{x)yo{x), in (0,1). 

Using once more the Theorem[5]-ii) and taking into account that v is nonpositive 
and yo is nonnegative in [0, 1], we get 

yt{x,t) < max (yl^{x) + v{x)yoix))+ < max (2/^'(a;))+ V {x,t) e Qt- 
xe[o,i] xe[o,i] 

The lower bound of b) is a direct consequence of the definition, the homoge- 
neous boundary condition and that y is nonnegative. To derive the upper bound, 
we will use Bernstein's method (see [14j p. 414 and 537]), by introducing the 
auxiliary function w{x,t) = e^^"^'*) + pe^~^ — 1 with p = max2,g[o,i](j/o(2;)e^°*^'^^). 
Let us remark that p > 0, thanks to the hypotheses on yo- 

By direct calculation it can checked that w G C^'^{Qt) and satisfies the 
following PDE: 

wt{x,t)-wxxix,t) = {v{x)y{x,t) - iyxix,t)f) ey'^^''^ - pe'-\ 
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Let us denote by {xo,to) a point in Qt where w takes its maximum value. We 
will see that xq — 0. If {xo,to) G Qt, it is well known that Wt{xQ,to) — and 
Wxxixoyto) < 0, and using the previous PDE we arrive to the contradiction 

< Wt{xo,to) - w^^{xo,to) < -pe^"^« < 0. 

The same argumentation can be used to exclude the case to = T, because the 
unique difference is that Wt(a;o,io) > 0. Consequently, {x(),to) must be a point 
of the "parabolic" boundary of Qt, that is, = 0, = 1 or to = 0. The last 
two possibilities can be excluded by noticing that w{0,t) = pe > p — 
w{x, 0) = e^"^^^ + pe^~^ — 1, and taking into account the choice for p: 

w^x, 0) = y[,(a;)e^''(^) - pe^"" <p- pe^^^ < 0, Vx G [0, 1]. 

Therefore, 0) is non-increasing in [0, 1] and w{x,0) < w(0,0) — pe, for all 
X G [0,1]. Consequently, w{x,t) < pe — 'w{0,t), for all {x,t) £ Qt and hence, 
Wx{0,t) < for all t £ [0, 1]. This is equivalent to say that yx{0,t) < pe for all 

[0,1]. 

The proof for c) is similar, using now the auxiliary function 'w{x, t) = e''^^'*) — 
pe^ — 1 with 

p= min (y(|(a;)e2'«(^)) < O.B 

KG [0,1] 

COROLLARY 1 Under the same hypotheses of Proposition QJ except that 
V £ L°°{0, 1), the conclusions a) — c) remain valid, for almost every point. 

Proof of Corollary [TJ By the classical convolution technique, it can be 
deduced the existence of a sequence {wfcl^i C C[0, 1] such that — ||w||ioo(o,i) < 
Vk{x) < for all fc = 1, 2, . . . , in [0, 1] and Vk — ?> w in L'(0, 1), as fc — > +00, for 
some q > 2. If we denote by y^^ £ C^'^{Qt) the classical solution of problem 
(|4.28l) with V — Vk, applying 14, Theorem 9.1, p. 341] to the difference y^,^ — y 
and using that < yv^,{x,t) < || J/oI1l°°(o,i) in Qt, by Theorem[5l we get 

\\yv,-y\\w^-\QT)-^\\'"''~'"\\L''(QT)~^^ A:^-+oo. 
In particular, 

iyvk)t — >yt vl\L''{Qt), as fc-j'+cxD. 

Combining the estimate above with [Ml Lemma 3.4, p. 82] for the values m = 
s = 1, r = and q > 2, we deduce 

(y«J:r(0,t) 2/^(0, t) and (2/^J^(l,i) ^ y^(l,t) in L«(0,r), as fc ^ +00. 

Using now the corresponding expressions a) — c) for each y^^. and the fact that 
the lower and upper bounds do not depend on Vk, we can pass to the limit in 
them and arrive to the same ones for y that will be satisfied for almost every 
point. ■ 
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4.2.4 Proof of Theorem [3 

The conclusion of the Theorem [7] is clearly valid if the L^— norm of the initial 
datum is small enough compared with e. More precisely, if |jyollL2(o.i) ^ f '^"^'^ j 
it is sufficient to take Ti = T2 — . . . ~ Tm — and mi — m2 = . . . = jum = 0. 
Of course, we want to deal with the general case. 

Let us begin with a useful observation. As we have mentioned before, it 
is well known that for each j/o G i+(0,l), we can construct an approximat- 
ing sequence {yok}kLi C C^+°'[0, 1] with compact support in (0,1) such that 
yok{x) > for all x G [0, 1] and k, with yofe — >■ yo in L^{0, 1), as fc — >■ +00. If we 
denote by yy^ the unique solution of (I4.28|) and by yy^i^ the unique solution of 
the same problem with initial datum j/ofc instead of yo, subtracting both PDE, 
multiplying it by the difference yyg — yy^,, and integrating by parts, we arrive to 
the expression 

^ {yyo(.x,T)-yy„,{x,T)fdx~^J^ (jjo^x) ~ yok{x)f dx = 



iiyvo(.x^t)-yyokix,t))x) dxdt+ v{x){yy„{x,t)-yy„^{x,t)) dxdt<0, 



Qt 

thanks to the hypothesis v{x) < 0, a.e. x e (0, 1). Therefore, 

1 .1 
{yvo{x,T) - yy^i^{x,T))'^dx < / {yo{x) - yok{x))'^dx, 
Jo 

and, consequently, it is possible to make the L^—norm of the difference of the 
solutions at time T as small as needed, by taking the L^— norm of the initial 
conditions small enough. As a consequence of this observation and taking into 
account the conclusion of the Theorem [7] that we want to derive, we will assume 
in the sequel that the initial data of the problems appearing along the proof are 
nonnegative regular functions with compact support in (0, 1). 

Now, we are in conditions to prove the existence of the nonnegative constants 
TOi, . . . ,tomj the intermediate times < Ti < r2 < . . . < Ti\i < T and the 
corresponding solutions yj for j = 1, . . . , Af (see (|4.17p - (|4.19p ). such that the 
conclusion of the Theorem [7] holds. 

For that purpose, we consider the family of problems 

yt{x,t) = y^^{x,t) ~my{x,t)x{o,i){x), in Qt, 
y{0,t) = 2/(l,t) = 0, in (0,r), (4.30) 

y{x,o) = yoix), in (0,1), 

and denote its solution by ym,i 

for each m > 0. Multiplying the PDE by ym,i 
and integrating by parts in Qt, we get as above 



{ym,i{x,T)) dx~- / {yo{x)) dx 



■10 ^ Jo 
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Consequently, 



{{ym,i)x{x,t)fdxdt-'m / {ym,i{x,t))'^dxdt. 
Jo Jo 



{ymAx,t)Ydxdt<— I {yo{x)ydx. (4.31) 



From this inequality we deduce that the sequence gm{t) = y jlj^{ym^i{x,t)ydx 
converges towards as m — ^ +00 in L^(0,T) and (taking a subsequence, if 
necessary) we have that gm{t) ^ as m — > +00 for a.e. t g (0,T). So, given 
e > 0, there exists Ti e (0,T) and toi > such that ^^^(ri) < j-^==. 
Denoting j/i — ymi.i this is equivalent to say 

'(yi(-,TO)^d.<^^^^. (4.32) 

Following the same scheme, we consider now the family of problems 

yt{x,t) = y^^{x,t) - ■my{x,t)x{i,2i){x), in (0, 1) x (Ti, T), 
2/(0, t) = 2/(l,t)=0, in (Ti,r), (4.33) 

y{x,Ti) = yoi{x), in (0,1), 

with yoi{x) = yi{x,Ti), yi being the function obtained previously. If we denote 
its solution as ym,2 for each m > 0, repeating exactly the same argumentation, 
it is possible to derive the existence of T2 € (Ti,T) and m2 ^ such that 
J/2 = J/m2,2 satisfies 

21 2 
(^^(^'^^))'^^^4(2]^- (^-3^) 
In fact, we can select T2 as close to Ti as needed: more precisely, satisfying 

0<r2-T,<min{ ^^^/^^^^ ,T-T,} (4.35) 

with 

Ci = ||yoi||L~(04) (l ■ ma^(y'^,{x))+ + e • max (y^,i(a;)e^-(-))) . (4.36) 

At this point it is evident that we are selecting the initial datum yoi to be a C^"'"°' 
function, with a S (0, 1), with compact support in (0, 1) (see the beginning of 
this section). This is needed for the right-hand term of (14. 36^ to be well defined. 
Let us remark that the terms appearing in (j4.36p are clearly related with the 
conclusions of Proposition [1] (items a) — b)). 

Under these conditions, we want to prove that 

"(y2(x,T2))^dx<^^^|^. (4.37) 
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Taking into account (|4.34l) . this will be true if we derive 

'{y.M)fd.<^^^^^. (4.38) 

This is the most dehcate point of the proof, where Corollary[T]is essentiaL Again, 
we multiply the PDE of problem ()4.33|) with m = m2 by 1/2 and integrate by 
parts in the domain {x,t) e (0,1) x {Ti,T2), to obtain 

^ {y2{x,T2)fdx-~ {yoi{x)fdx = 
/ {{y2)^{x,t)ydxdt+ / {y2)^{l,t)y2{l,t)dt, (4.39) 

/Ti Jo JTi 

where the boundary condition at a; = has been utilized together with the fact 
that (0, Z) n (1,21) — 0. In particular, combining this relation with ()4.32p . we 
have 

iy2ix,T2)fdx < [ {yoiix)fdx + 2 [ {y2)x{l.t)y2{l.t)dt < 



Ti 



T2 



and the estimate (|4.38p will hold, if we are able to prove that 

^^y.).(M)..(M)rfi<^(^^. (4.41) 

Let us show that this holds thanks to the choice of T2 (see (|4.35p - (|4.36p ). Using 
Corollary [1] (with v{x) ~ —ni2 ■ Xii,2i)ix) that is only a bounded function) for 
problem (|4.33p . it follows that 



(2/2)x(^0 = (2/2)^(0, 0+ / {y2)xxix,t)dx^iy2)A0,t)+ / {y2)t{x,t)dx < 

Jo Jq 

< e • max (y' (x)e^°i(^)) + 1- max (2/01 (x))+ C2. 

We can now combine this estimate with the fact that < y2 l£ l|yoi||L°°(o,i) 
thanks to the Theorem [S]-i) and ii), to get 

T2 /.T2 2 

^ - 4(2A/ - 1) ' 



(.y2)Al,t)y2{l,t)dt < 2 / C2-\\yai\\L^io.i)dt = 2{T2-T{)C. < 

that is exactly (|4.4ip . 

Finally, we can repeat the argumentation to select the constant rrij , the time 
Tj e [Tj-i, T) and the solutions yj, j — 1,2, ... ,M — 1, verifying 
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|(,,(.,T,))^d.<||^, (4.42) 

and arrive to 

/■I ,2 

J^{yM{x,TM)fdx<-, (4.43) 

as desired. ■ 
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